Hierarchy of relevant couplings in perturbative renormalization group transformations 
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The phase diagram for the interacting fermions in weak coupling is described by the perturbative 
renormalization group equations. Due to the lack of analytic solutions for these coupled non-linear 
differential equations, it is rather subtle to tell which couplings are relevant or irrelevant. We 
propose a powerful classification scheme to build up the hierarchy of the relevant couplings by a 
scaling Ansatz found numerically. To demonstrate its superiority over the conventional classification 
for the relevant couplings, we apply this scheme to a controversial phase transition in the two-leg 
ladder and show that it should be a non-trivial crossover instead. The scaling Ansatz we propose 
here can classify the relevant couplings in hierarchical order without any ambiguity and can improve 
significantly how we interpret the numerical outcomes in general renormalization group methods. 
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Renormalization group [3, 0, 0, B| is a powerful method 
to determine the effective interactions for the complex 
system in low-energy limit. By integrating out the de- 
grees of freedom at the longer length scale, the couplings 
describing the effective interactions flow according to a 
set of RG equations. This approach has been success- 
fully applied to a wide variety of physical phenomena, 
including the transport theory in the presence of im- 
purity scattering 1], ground-state properties for cold 
atoms [3 B , phase diagrams for iron pnictidesQ, Kondo 
lattice[10[ and other correlated systems. RG is particu- 
larly helpful when the quantum fluctuations in the sys- 
tem are strong 
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so that the 

mean-field description is invalid or hard to justify. For 
instance, it has been demonstrated that, despite of the 
repulsive interactions at the short length scale, electrons 
form pairs with unconventional rf-wave symmetry in the 
low-energy limit in the two- leg ladder [l3. Ilij . In addi- 
tion, the RG analysis also predicts that the carbon nan- 
otubes are Mott insulating spin liquid [H, S3, EH ElEl , 
which is beautifully realized in recent experiments [2J]. 

In weak coupling, the phase diagram for the ground 
states of the interacting fermions is described by the cou- 
pled non-linear RG equations, 
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where §i are the couplings for the effective interactions 
and I is the logarithmic length scale. Ay denotes the 
scaling dimensions at the tree level and Bj kl are 

the one-loop and the two-loop renormalization. While 
the derivation of the RG equations are rather standard 
and can be found in the literature 0, E3, E3], inter- 
preting the RG flows obtained in numerics can be sub- 
tle and tricky. When the tree-level contributions are 
not zero, reading off the relevant couplings are rather 
straightforward [271 \2ty . The challenge arises when all 
couplings become marginal, i.e. Ay = 0. The standard 
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FIG. 1: RG flows of the eight couplings in the doped two-leg 
ladder. The initial couplings are generated from the on-site 
interaction U/t = 10 -5 and the velocity ratio V2/V1 = 7.8. 
We only plot the flows within the perturbative regime where 
all couplings gi < 1. 



recipe is to integrate the RG equations up to the cutoff 
length scale l c where the maximal coupling is of order 
one. At I — l c , we identify the couplings gi(l c ) ~ 
to be relevant while those couplings gi(l c ) C 1 as irrel- 
evant. But, this face-value classification scheme at the 
cutoff length scale does not always work and some cou- 
plings are ambiguous to be identified as either relevant 
or irrelevant. 

To amend this ambiguity, Ledermann, Le Hur, and 
Rice[29j come up with a clever method to patch the RG 
flows at different length scales and obtain a hierarchy of 
gap opening due to different Fermi velocities. Inspired 
by their success in building up the hierarchy of relevance 
in the RG flows, we analyze the RG equations for corre- 
lated electron systems in weak coupling numerically. Af- 
ter extensive numerical studies, we find that the relevant 
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FIG. 2: Exponents extracted from the numerical solutions for the RG equations at different velocity ratios V2/V1 in the scaling 
regime. For small velocity ratios, all couplings except /f 2 are equally dominant with exponents ji — 1. As the velocity ratio 
increases, the relevance of the couplings can be classified by the exponents accordingly. 



couplings can be well captured by the scaling Ansatz 14 1 
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revisit the phase diag ram of the doped two-leg ladder. 
In previous study [l2j, it was shown that the different 
Fermi velocities can lead to a quantum phase transition 
from the spin liquid to another gapless phase. How- 



where Gi are some non-universal constants and Id is the 
divergent length scale from the one-loop RG equations. 
The relevant couplings are thus captured by the expo- 
nents < 7^ < 1 appearing in the scaling Ansatz as 
shown in Fig. 1. It may sound odd at first glance why 
the divergent length scale Id appears in the scaling Ansatz 
for RG flows in the perturbative regime. This is due to 
an approximate scaling relation in weak coupling which 
we will come back later. 

To demonstrate the validity of the scaling Ansatz, we 
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ever, in later numerical studies[30|, [31|, there is no hint 
for the phase transition. Therefore, it remains contro- 
versial whether the difference in Fermi velocities can 
drive a quantum phase transition. It was shown be- 
fore that there are eight independent couplings, g = 

( C ll 5 c i 
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/i2> fi2)i describing the forward 
and Cooper scattering in the system. Readers who are 
interested in details are encouraged to read the previous 



publications 12, 2 



After coarse-graining the fluctuations at shorter length scale, the one-loop RG equations for the couplings are 

k 
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where Vi are the Fermi velocities and the short-hand notations arc defined, hij = 2ftj + 6ijCu and ciijk = (vi + Vk){vk + 
Vj)/2vk(Vi + Vj). ^ 



A typical example for the RG flows with the on-site in- 
teraction U/t = 10~ 5 and the velocity ratio V2/V1 = 7.8 
is shown in Fig. [T] First of all, we can extract the diver- 
gent length scale Id ~ 10 7 numerically. At the beginning, 
the RG flows are rather messy and it is hard to tell which 
are the relevant couplings. However, as the short-range 



fluctuations are integrated out gradually, the couplings 
enter the scaling regime described by Eq.([2]) with differ- 
ent exponents 7,. The straight lines in Fig. [1] indicate 
that the scaling Ansatz works extremely well. It is worth 
emphasizing that all renormalized couplings remain in 
the perturbative regime even though the Ansatz myste- 
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riously contains the divergent length scale Id- 

The exponents provide an alternative way to classify 
the relevant couplings. The eight couplings fall into four 
categories: (1) the predominant group 7^ = 1, (2) the 
first subdominant group 7 c f> 2 = 7^ =3/4, (3) the second 
subdominant group 7^ = 7 cfi = -y^ = 7^ = 1/2, 
(4) the subsidiary group jf* = 0. By extracting the 
exponents from the trends of the RG flows, we can build 
up the hierarchy of relevance without any ambiguity. 

The classification scheme can be applied to different ve- 
locity ratios with exponents summarized in Fig. [5] Upon 
doping away from the half filling, the velocity ratio Vijv\ 
starts from unity and gradually increases. When the ve- 
locity ratio is close to one, all exponents 7^ = 1 except 
that for fi 2 - It is not surprising that the conventional 
face- value classification scheme works rather well in this 
regime since all couplings are equally relevant. Apply- 
ing the bosonization technique, the interactions related 
to gap opening are 
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Here we express the Dirac fermion fields in term of 



chiral boson fields, ippia = Ki a e 



with the Klein 



factors k,vv[16|. Meanwhile, by introducing the conju- 
gate boson fields for each flavor, Lp ia = 4>Bia + 't'Lia 
and 9i a — 4>m a — 4>Lia, we separate the boson field 
into charge and spin modes, 9 ip = (0^ +%) /V2 and 
®i<j = (#it — /v2 (similar for the conjugate field ip). 
Since the couplings c^, C22, cj 2 are relevant, spin gaps Ai 
and A2 for the antibonding and bonding bands develop 
under RG transformation. The magnitudes of the spin 
gaps can be estimated by the semiclassical approxima- 
tion, Ai/A 2 ~ flcfj + c\ 2 ) I (|cf 2 | + C12), where the val- 
ues of the couplings are taken at the cutoff length scale. 
The detail bosonization analysis shows that the ground 
state is a spin liquid with unconventional <i-wave pairing 
between electrons 12, TBI. Hij. 

As the velocity ratio increases i^/fi > 5, 7 c j 2 remains 
one but the other exponents start to fall, as shown in 
Fig. [21 The four-level hierarchy of relevance described 
before starts to emerge. When the velocity ratio further 
increases vilv\ > 8, the second subdominant exponents 
also falls to zero, 7 C <> = 7^ = 7 c p 2 = 7^ = but 
the predominant and the first subdominant exponents 
remain intact. The exponents of the scaling Ansatz pro- 
vide a simple and elegant way to classify the relevant 
couplings. It is also interesting that the evolution of the 
exponents mainly stays on plateaus of simple ratio num- 
bers. These ratio numbers can be obtained by performing 
linear stability analysis [26( near the fixed rays of the RG 
flows. However, since the initial couplings are nowhere 
close to the fixed rays, it is not clear why the RG flows 
inherit the simple rational exponents. The plateau struc- 
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FIG. 3: RG flows up to two- loop order for the charge sector 
in the one dimensional Hubbard model. 



ture of the exponents seen in numerics remain open for 
further investigations. 

Within our classification scheme, we find the pattern 
for the excitation gaps is the same except the magnitude 
of the spin gap in the bonding band A2 is much larger 
than Ai in the antibonding band. Therefore, as the ve- 
locity ratio 1)2 /vi increases, the gap ratio A2/A1 also 
increases significantly but there is no phase transition 
for 1 < V2/V1 < 8.6. However, the conventional read-off 
scheme at the cutoff length scale misses all the subdom- 
inant couplings and misidentify the crossover as a quan- 
tum phase transitionflil. Il5j |. In addition, the scheme 
proposed here naturally leads to a hierarchy of excita- 
tion gaps predicted by the patched RG method [2^]. 

We also try out the scaling Ansatz beyond the one-loop 
order and found it remains valid. For instance, the RG 
equations for the one dimensional Hubbard model[32j are 
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Here we only put down the charge sector since the RG 
equations for the spin sector are decoupled with similar 
structure. As is evident from Fig. [31 the two-loop cor- 
rections are negligible and both G p and G c follow the 
scaling Ansatz rather well. Therefore, the scaling Ansatz 
and the corresponding exponents remain valid beyond 
the one-loop order as long as the flows stay in the per- 
turbative regime. 

Now we come back to the puzzle about the diver- 
gent length scale Id in the scaling Ansatz. This arises 
from a non-trivial relation for the one- loop RG flows, 
gi{l) = UGi(Ul), where Gi(l) is the solution with order- 
one initial values while gt{l) is the solution in weak cou- 
pling (with initial values of order d « 1). Thus, there 



4 



is a non-trivial connection between the RG flows in per- 
turbative regime and those in the singular regime. This 
relation may seem an artifact for the one-loop RG equa- 
tions. But, it is rather remarkable that the relation re- 
main approximately correct even when higher-loop cor- 
rections are included. This is also consistent with the 
numerical findings that the RG flows are dominated by 
the one-loop terms with negligible corrections from the 
higher loops. 

Now we try to see how the scaling Ansatz emerges in 
the non-linear RG equations. It is insightful to rewrite 
the equations in matrix form dgi/dl — g T AiQ 1 where 
the real symmetric matrices are (Ai)jk — A\ . Assume 
that the matrices Ai, i — 1,2, ..,n, are positive definite. 
One can use the standard decomposition to rewrite the 
matrix A n = L n D n L^, where L n is a unit n X n lower- 
triangular matrix and D n is a diagonal Gaussian pivot. 
Simple algebra leads to the important inequality, 

where A„ > is the n-th pivot in the diagonal matrix 
D n . Integrating the inequality in Eq. ([5]), one can show 
that the RG flows become singular at some finite length 
scale I — Id- Following similar algebra, one can also show 
that all couplings become singular at exactly the same 
length scale. 

The above proof establish the existence of the singular 
length scale Id in mathematical rigor, making the scal- 
ing Ansatz with power-law singularity plausible but not 
yet proven. However, if one plugs in the scaling Ansatz 
into the inequality, it leads to the constraint ji < 1 as 
found in the numerics. The situation is further compli- 
cated by the fact that the matrices Ai in the RG equa- 
tions are not always positive definite. This is reasonable 
because, for some initial coupling profile, the RG flows 
may not become singular at all. However, as long as the 
initial couplings sit inside the unstable manifold in the 
multi-dimensional coupling space, the above statements 
are expected to be correct. 

Since the classification scheme we proposed here is 
rather general, the numerical solutions from the func- 
tional RG approach [2^ can be classified by the scaling 
exponents as well. In fact, we apply the scaling Ansatz 
to a recent RG analysis |33[ for iron pnictides. There are 
four couplings (u±, 112, U3, U4) describing the interactions 
between the electron and hole pockets. We numerically 
integrate these RG equations and find the scaling Ansatz 
works again with exponents ji = 1, - , 1, 1. 

In conclusion, we propose a powerful scheme to clas- 
sify all relevant couplings in the one-loop RG equations. 
As long as the initial couplings are weak, the renormal- 
ized couplings enter the scaling regime characterized by 
a unique set of exponents 7*. We can build up a hierar- 
chy of relevance by these exponents without ambiguity 



and greatly improve the interpretations of the results ob- 
tained by various RG methods. 
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